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ABSTRACT 

Dynamical friction arises from the interaction of a perturber and the gravitational 
wake it excites in the ambient medium. This interaction is usually derived assuming 
that the perturber has a constant velocity. In realistic situations, motion is accelerated 
as for instance by dynamical friction itself. Here, we study the effect of acceleration 
on the dynamical friction force. We characterize the density enhancement associated 
with a constantly accelerating perturber with rectilinear motion in an infinite homo- 
geneous gaseous medium and show that dynamical friction is not a local force and 
that its amplitude may depend on the perturber's initial velocity. The force on an ac- 
celerating perturber is maximal between Mach 1 and Mach 2, where it is smaller than 
the corresponding uniform motion friction. In the limit where the perturber's size is 
much smaller than the distance needed to change the Mach number by unity through 
acceleration, a subsonic perturber feels a force similar to uniform motion friction only 
if its past history does not include supersonic episodes. Once an accelerating perturber 
reaches large supersonic speeds, accelerated motion friction is marginally stronger than 
uniform motion friction. The force on a decelerating supersonic perturber is weaker 
than uniform motion friction as the velocity decreases to a few times the sound speed. 
Dynamical friction on a decelerating subsonic perturber with an initial Mach number 
larger than 2 is much larger than uniform motion friction and tends to a finite value 
as the velocity vanishes in contrast to uniform motion friction. 

Key words: hydrodynamics ~ ISM: general 



1 INTRODUCTION 

Dynamical friction is the reaction force that a perturber feels from the gravitational wake it generates in the ambient medium. 
This fundamental interaction has therefore had a wide range of applications in the formation of astronomical structures. 
Applications to coUisionless systems include mass segregation in star clusters (e.g. Kim and Morris 2003, Portegies et al 
2004), galaxy distribution in hierarchical cluster formation (e.g. Frenk et al 1996, Bullock and Johnston 2005), the decay of 
satellites in galactic halos and galaxy mergers (Binney and Tremaine 2008 and references therein). Dynamical friction also 
occurs in collisional particle systems and was applied to the study of the growth of planetesimals (Hornung et al 1985, Stewart 
and Whetherill 1988), the eccentricity excitation of planetary embryo orbits (Ida 1990, Kokubo and Ida 1996, Namouni et 
al 1996), and the confinement of planetary rings (Goldreich and Tremaine 1982). For gaseous systems, applications include 
the orbital decay of compact stars orbiting supermassive black holes (Narayan 2000), and the heating of intercluster gas by 
decelerating supersonic galaxies (El-Zant et al. 2004, Faltenbacher et al 2005, Kim et al 2005). 

Following the pioneering work of Chandrasekhar (1943) on the gravitational drag experienced by a perturber in coUisionless 
star systems, the derivations of dynamical friction assume that the perturber has a constant velocity in the ambient medium. 
The assumption of uniform motion is made for coUisionless systems (Binney and Tremaine 2008), collisional particle systems 
(Hornung et al 1985, Stewart and Whetherill 1988, Ida 1990, Goldreich and Tremaine 1982) as well as gaseous media (Kim 
and Kim 2007, Kim et al 2008, Ostriker 1999). The resulting dynamical friction force which depends on the perturber's 



* E-mail: namouni@obs-nice.fr 



2 Fathi Namouni 



velocity is subsequently used to model the evolution of the perturber's motion whose velocity varies with time. Only the 
numerical simulations of the nonlinear fluid equations are free from this assumption (Sanchez-Salcedo and Brandenburg 1999, 
2001). The neglect of the perturber's acceleration in dynamical friction is justified by the small amplitude of the force that is 
proportional to the square of the pertuber's mass. Neglecting acceleration is independent of the type of motion and is mainly 
motivated by the lack of an analytical framework where the perturber's accelerated motion could be taken into account. To 
overcome this obstacle, we derive in Appendix A a new expression of dynamical friction that applies to a general type of 
motion and depends on the perturber's motion and the boundary conditions but not the density perturbation. In particular, 
dynamical friction is shown to be non-local as it depends on the whole trajectory of the perturber. Although non-locality 
does not appear if the perturber's velocity is constant, its origin is easy to identify: as the perturber's wake propagates in 
the medium, a velocity change because of acceleration affects the density enhancement near the perturber and offsets the 
force from its local value as the far parts of the wake were generated by a different velocity. Therefore when acceleration is 
taken into account, the amplitude of dynamical friction may change with respect to the value derived assuming a constant 
velocity because in modelling the force on an accelerated perturber as that on a constant velocity perturber, one assumes that 
the gravitational wake produced by accelerated motion is similar to that of uniform motion. We shall see that this difference 
between the two wakes makes dynamical friction dependent on the initial velocity of the perturber. Another difference between 
the wakes generated by constant velocity motion and accelerated motion is the absence of a perturbed region with front-back 
symmetry with respect to a subsonic accelerating perturber. In effect a constant velocity perturber raises a gravitational 
wake in the ambient medium that is symmetric with respect to it and truncated by the sonic Shockwave. This truncation 
defines two regions, one that lacks front-back symmetry located behind the perturber, and one with front-back symmetry 
around the perturber that becomes the largest region inside the wavefront for velocities smaller than half the sound speed. 
Dynamical friction arises from the former region as shown by Ostriker (1999). When acceleration is taken into account, front- 
back symmetry is lost everywhere. This would intuitively imply that dynamical friction on an accelerated perturber should 
be much larger than the value derived assuming uniform motion for subsonic perturbers. Curiously, we shall show that the 
force on an accelerating subsonic perturber initially at rest is similar to the uniform friction force. Lastly, accelerated motion 
friction would differ from uniform motion friction because the former never reaches steady state. In effect, dynamical friction 
on a constant subsonic velocity perturber is known to be time- independent (Ostriker 1999). 

It is the object of this paper to quantify how taking into account acceleration modifies dynamical friction on a perturber 
travelling in an infinite homogeneous gaseous medium. The perturber's motion is assumed to have a constant acceleration, a 
choice that makes the analytical developments simpler. Physically such constant acceleration approximates the deceleration 
felt by the perturber because of dynamical friction as the force amplitude is small and the corresponding instantaneous velocity 
change is equally small. Constant acceleration also models the motion of a young star with an asymmetric jet system moving in 
interstellar gas (Hirth et al 1994, Hartigan and Morse 2007, Namouni 2005, 2007). In what follows, we characterize the density 
perturbation associated with accelerated motion in Section 2. We derive the force on an accelerated perturber in Section 3 
using the general expression of dynamical friction (derived in Appendix A) that does not depend on the density enhancement 
caused by the perturber but only on the perturber's motion and the medium's boundary conditions. In particular we show 
that including acceleration makes dynamical friction dependent on the initial velocity of the perturber. Section 4 contains 
concluding remarks. In Appendix B, the new expression of dynamical friction derived in Appendix A is validated by applying 
it to uniform rectilinear motion and recovering the known results on dynamical friction. 



2 DENSITY ENHANCEMENT 

Consider a point-like perturber of mass M on an accelerated rectilinear trajectory, of initial velocity Vb and acceleration A, 
that travels inside an inviscid isothermal gaseous medium of density p and sound speed c. Assuming the density perturbation 
associated with the perturber's wake is small compared to the medium's unperturbed density po, the corresponding density 
enhancement q defined as p = po(l + f?) satisfies the standard forced sound wave equation (Appendix A, equation IA7|l : 

V'q - ^ dh = -^VVp = -47r^ 5'[x mWit), (1) 

where </>p — ~GM/\x — i{t)\ is the perturber's gravitational potential, J^{t) is the Heaviside function, G is the gravitational 
constant, ^(t) = {At'^ /2 + Vot)ez is the perturber's trajectory chosen along the z-direction and 5{x) is Dirac's delta function. 
The solution of this equation is the Lienard-Wiechert potential which may be written in integral form using the retarded 
Green function as 



The acceleration A introduces two natural time and length scales given by: 
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T-^,L^2_. (3) 

In these units, the perturber's trajectory may be written as ^(t) = + ./#oi where ./#o is the initial Mach number. The 
Mach number is a function of time and reads: ^ = t + ^^q. The scaling time T is the duration to increase the Mach number 
by unity and the scaling length L is the corresponding travelled distance. Applied to the density enhancement these scalings 
yield: 



5 ( u-t+ln^ + (z- .^ou-u^/2YY^] .je{u) 

. . GMA I- y [ y " ' n J ' 

g{x,t) = —- d^, (4) 

C 7_oo + (2 - ^OU - U2 /2f] ' 



where B? = x^+y^ is the cylindrical radius; unless it is stated otherwise, all length and time variables and constants are normal- 
ized to (l3|. In order to derive the density profile we use the property of the 5- function a(u)5[/3(u)] = '^^a{ui)5{u — Ui) /\fi' {ui)\ 
where Ui are the roots of Piu), and /3' is the derivative of /3. The roots of the argument of the 5-function in equation Q 
satisfy the quartic equation: 

+ ^ou^ ~{l + z~ J(l) u + 2(t - Mqz)u - + r-^ = 0, (5) 

where r — {R^ + z^Y^'^ is the spherical radius. This equation is the counterpart of the quadratic equation derived by Ostriker 
(1999) for uniform motion. The former may be reduced to the latter once time and length dimensions are restored. Equation 
([5|, however, cannot be solved analytically because its varying coefficients do not permit simple and unique expressions of 
the roots. Instead we determine the roots numerically and substitute them into the density expression obtained by applying 
the 5-function property cited above. The density expression reduces to: 

Q{x,t) = \t - - {.^0 + Ui){z - .^QUi - /2)\~'^. (6) 



roots u 



of cq. j5t 



Fig. (1) follows the level curves of g for perturbers initially at rest (./#o = 0) as they evolve from t = to Mach 6. The front-back 
symmetry that characterizes the density enhancement of subsonic constant velocity motion in a uniform medium is absent 
for both subsonic and supersonic motion. The perturber reaches Mach 1 at z = 0.5 but remains inside the sonic wavefront of 
radius t as the latter has been propagating faster than the perturber's velocity. Once the velocity is supersonic, a Mach cone 
is formed inside the sonic wavefront and a shallow density depression develops and trails the perturber {z ^ —1). Eventually 
as Mach 2 is reached the perturber catches up with the sonic Shockwave and exits the sonic region to form an external Mach 
cone. Fig. (2) shows the level curves of g associated with decelerating perturbers with 7^ as a function time until the 
perturber's velocity vanishes. The density enhancement of subsonic perturbers has a structure similar to that of uniform 
motion (Ostriker 1999) but lacks its front-back symmetry with respect to the perturber. We show in the following section 
that under certain circumstances, dynamical friction may be not significantly affected by the loss of front-back symmetry even 
as closed level curves around the perturber contribute to the force unlike those of uniform motion (see Introduction). The 
deceleration of supersonic perturbers modifies the structure of the density enhancement and reduces the extent of the Mach 
cone but does not erase it as the velocity vanishes (bottom row, right panel) indicating the presence of a finite friction force. 



3 FRICTION FORCE 

The reaction force exerted on a point-like perturber in a homogeneous gaseous medium by its gravitational wake is given by 
equation (IA9|) and is written in dimensional variables as: 

^ J aVly+i(t-r/c)],r^ct If "I 

where r, 6 and ip are the spherical coordinates of the vector y defined as the relative position within the medium with 
respect to the perturber's position at the retarded time t — r/c. It may be written as y = x — ^{t — r /c) where x is the 
position of fluid element in the medium. This definition of y influences the force integral through the boundary conditions 
of the medium denoted by dV. The phase A — ^(t) — ^{t — r/c), the upper radial cutoff at the retarded sonic Shockwave 
r ^ ct, and the boundary conditions are all that is necessary to compute the friction force. In particular, the absence of an 
explicit expression for the density enhancement is irrelevant for the derivation of the force. The derivation of expression ((TJ 
assumes that the perturber is assembled at t = a fact that introduces the size of the retarded sonic Shockwave as an upper 
cutoff. This assumption is not restrictive as it allows us to study the time- dependent evolution of dynamical friction as well 
approximate steady state behaviour (see section 3.2). It also describes turning on a numerical simulation of dynamical friction 
(Sanchez-Salcedo and Brandenburg 1999, 2001). 
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For rectilinear motion along the ^-direction in an infinite medium, the force is equally along and its expression can 
be reduced further by integrating over the longitude and the variable r = cos 6 as: 

2'Kje{t){GMfpo f dr r r - r A 



F 



aV[y+i{t-r)],,%t 



dr 
A2 



1/2 



(8) 



ll + (A/r)2-2(A/r)r| 

where the scalings ([3]) have been applied; the perturber's trajectory ^{t), the time, t, and the position x and its spherical 
components are dimensionless hereafter. The expressions of Tmin and Tmax as a function of the perturber's trajectory and 
the variable r are determined from the boundary conditions. For uniform rectilinear motion, ^ is constant and the phase 
A/r = Calling ro the gravitational potential's small distance cutoff scaled to L ((3]), the integral that appears in the 
expression ((8]l may be written as (Appendix B, equation IB6I Ostriker 1999): 

(1 



lo = 



[1 + sign(l - .-#)] - log 



1 + . 



- [1 - sign(l - .^)] log 



ro 



(9) 



For accelerated motion, the phase A/r — + t~ r/2. In order to examine how acceleration modifies the dynamical friction 
force, we consider the three cases of perturbers initially at rest {^o = 0), accelerating perturbers (./#o > 0) and decelerating 
perturbers {.y£o < 0). 

Before deriving dynamical friction for accelerated motion, a choice of the small distance cutoff ro must be made. This 
cutoff is usually taken as the perturber's accretion radius or the system's typical size. The former dimensional cutoff reads 
2GM/V^ and is a function of the Mach number whereas the latter is a constant that does not depend on the fiow's variables. 
In this work, the perturber is treated as a compact object whose gravitational potential is truncated at small distances by 
the perturber's size. To satisfy such a condition and be consistent with a point-like potential, the perturber's size needs to be 
much smaller than the accretion radius 2GM/V^. This choice of cutoff is in agreement with the numerical simulation of the 
fiuid equations in the context of Bondi-Hoyle-Lyttleton accretion (Edgar 2005). For larger size or extended objects, equation 
((7| no longer applies and dynamical friction may be obtained from the general expression derived for an arbitrary perturbing 
potential (Appendix A, equation [ 



3.1 Perturbers initially at rest 

In an infinite medium, the only relevant physical boundary condition is that of the small distance cutoff \x — £,{t) \ ^ ro0 For 
a perturber initially at rest, £^{t) — 1? 12, lS./r — t — r/2 and the Mach number is ^ = t. The variables r and r therefore need 
to satisfy the three conditions: (i) -1 r ^ 1, (ii) r ^ [1 + (2f - r)^/4 - {ro/rf]/{2t - r) = f{r,t,ro), and (iii) r ^ t. The 
last condition implies that A > 0. Examination of the three conditions shows Tmin = —1, and that Tmax and the remaining 
bounds on r in the integral (|8]) depend on the value of t as the velocity increases. There are five time intervals that need to 
be considered: 



I 

II 
III 
IV 

V 



(1 + 2ro)^/' - 1< f 1 - (1 - 2ro)^/^ 
l-(l-2ro)^/2 <fs; l + (2ro)'/^ 
l + (2ro)^/^ ^ts; l + (l-2ro)'/^ 
1 + (1 - 2ro)'/' SC f 1 + (1 + 2ro)'''^ 



2ro) 



1/2 



sc t. 



Case I is that of transitory motion where the density perturbation is in the perturber's vicinity. Prior to t = (1 + 2ro)^''^ 
the force is zero as the three conditions above cannot be satisfied. The integration domain of r and Tmax are given as: 

r(^r^t with Tmax = /(r, t, ro), where r^ = 1 + t — [{1 + if — 2ro]^^^ . 

Calling / the integral in equation ([8]) and using these boundaries for r, we find: 

■r^ -2t'^r'^ + {4t^ - 4rot + 2r^)r 



2rlt 



2rot'^r{2t - r) 

ro(2 + r|) - (2 + m + r()t - rjr' + 
2rot2 



2t-ro 
2t^ 



log 



2t 



2t 



- ro 



2t3 



log 



2t 



(10) 

(11) 
(12) 



Examination of the variation of h as a function of time shows that it decays smoothly to zero as t = (1 + 2ro)^^^ — 1 is 
approached. The expression of the force in this time interval is the equivalent of that derived for uniform motion in Appendix 
B (equation IB7|) when ro/|l + ^| ^ ct ^ ro/|l — ^\ and the wake has not left the perturber's vicinity. In that case, the force 



^ The variables x, ^ and ro remain scaled to L ((Sjl. 
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read F oc ./#^^{— log[ct(l + ^/£)/ro] + ^ — ro/2ct + (1 — .^^)ct/2ro}; it helped explain the transient phase observed before 
saturation in the numerical simulation of the nonlinear fluid equations (Sanchez-Salcedo and Brandenburg 1999). 
Case II covers subsonic motion and the corresponding integration boundaries are: 

J^^^^'^ „,;<-v, - ,, = /(r,t,r-o), where r-J = t-l + 1)2 + 2ro]'''', 
.X = 1. 



,2 < r- < 
r1 



with 
with 



(13) 
(14) 



The latter domain (|14|) does not contribute to the force. The corresponding integrand is proportional to sign(l — t-\- r/2) — 1 
and the sign argument is positive for r > 2(t — 1). The expression of / for domain (|13|l may be obtained from equation (|lip 
by substituting rl for t. The result is ; 



In 



2n>t + (t - ro)(rJ + ri) ~ t^ijl - ri,) 2t - ro 



2rot2 



2t^ 



log 



rit + ro 



ro 



(15) 



In the limit ro <^ 1 and t < 1, Iu{ro <C 1) = t~^{2t — log |(1 +t)/{l — i)\). This is exactly the friction force acting on subsonic 
perturbers with uniform motion given by equation Q. It is interesting that although the density enhancements of subsonic 
accelerated motion and subsonic uniform motion differ significantly, the corresponding friction forces are identical when the 
limit ro <C 1 is taken for the accelerated motion friction. This result is the more surprising as the density enhancement 
associated with accelerated motion does not involve a region with front-back symmetry as explained in the Introduction 
leading one to guess wrongly that friction would be stronger for accelerated motion than it is for uniform motion. The 
explanation of this agreement comes from the friction force's dependence on the phase A = (_{t) — ^(t — r/c) (dimensional 
variables). As the contribution of the force for subsonic motion comes mostly from the vicinity of the perturber's retarded 
position {\y\ 0), we have r/c <C t and the phase becomes A = ^{t)r/c. The latter is exactly the phase of uniform motion 
with a Mach number ^ = C(^)/c regardless of the type of motion ^(t) as long as, (1) ro ^ 1 which implies in dimensional 
variables that |^(t)| <^ <? jro-, (2) the velocity remains subsonic and (3) the motion does not have a history with supersonic 
episodes (see section 3.3). 

Case III concerns motion between Mach 1 and Mach 2 where the perturber travels at supersonic speed but has not yet 
crossed the sonic shock wave. The integration domains are given as: 



r| ^ r < r2 with r^ax 

r^^r ^r% with 

ra ^ r < ri with 

ri ^ r ^ t with 



= f(r,t,ro), where r2 = t - 1 - [{t - l) - zvo 
= 1, where rl = t - I + [{t - if - 2ro\'^'^ 
= .f{r,t,ro), 
= 1. 



2rol^''2 



(16) 
(17) 
(18) 
(19) 



The integration radii satisfy r| < r2 r3 ^ 2{t — 1) ^ rl. This shows that domain (|19p does not contribute a finite term to 
the force in the same way as domain (I14|l . Domains (|16|l and p8[) contribute terms similar to that of equation (|lip except 
that for the former t is replaced with r2 and for the latter t is replaced with rj, and r| with r§. The second domain (|17p 
contributes the term J jp^ given as: 



(r - 2t)t f2 



log 



2t 



(20) 



which when combined with the contributions of domains (|16|l and psp yields 



/ill 



2rot - ro{r{ + r( - r^ + rg) + i(rj + r| + ri - r§) 2t + ro 



- + rj - r^ 
2ro 



2rot2 

2t-ro 
2t3 



2t3 



log 



rfi - ro 



■ ro 



log 



rft + ro 



ro 



(21) 



Case IV defines a small time interval near Mach 2 where the perturber exits the sonic wavefront. The corresponding integration 
domains are: 

r(^r^r2 with Tmax = /(r, ro), (22) 
r2 < r ^ r| with Tmax = 1, (23) 
rs < r < t with Tmax = /(r, t, ro). (24) 

The expression of / is obtained by summing up the terms associated with domains (|16p and (|17p of Case III, and substituting 



We remind the reader that the uniform motion friction force on subsonic perturbers Q is independent of ro and that no assumptions 
on ro are made to derive its expression. 
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t for rl in domain p8|) . The result is; 

, _ ro{rl-rl~ri-2) + {2 + ro + r9,-Tl + ri)t + {-r-i + Tl + ri)e -t' 



2t + ro 
2t3 



loe 



To 



2t 



ro 



2t3 



loE 



2t 



ro 



ro 



Case V concerns supersonic motion that exited the sonic wavefront > 2. The integration domains are given as: 



^2 ^r ^r2 
r2 ^r 



with Tmax = f{r,t,ro), 

with Tinax = 1- 



(25) 



(26) 
(27) 



The integral / is obtained as previously by combining the domains of Case IV and making the necessary boundary substitutions. 
The result is: 



ro{r2 



4) + {r^2 + rl)t~{r^2^ri)t^ 



2rot2 

For ro ^ 1 and t > 2, the scaled force reads: 



2t + ro 
2t3 



log 



2t + ro 



ro 



2t~ro 
2t3 



log 



2t-ro 



r{t- 



ro 



/v(ro < 1) = -r 



2(1 + log 2) + log 



(28) 



(29) 



This expression differs from that of equation ([S} through the first term in the bracket —2(1 + log2)/f^ implying that the force 
on accelerated motion of a perturber initially at rest, is larger than that on uniform motion. Fig. (3) shows the friction force 
(|12|) . (|15I) . (|2H) . (|25p . (|28|) as a function of time, or equivalently the Mach number, for the three values of the small distance 
cutoff ro = 10~^, 10~*, and 10~^. The force obtained by using the uniform motion expression, Jo, (O with ^ — t is also shown 
for the same values of ro. For subsonic motion, the agreement between the two forces is best when ro is smallest. The force Jo 
overestimates maximal drag by about 20%; it underestimates the drag on supersonic motion by the factor ~ 2(1 + log2)/f^. 
Lastly, we remark that the derivation of dynamical friction in the five time intervals shows that there are no force singularities 
associated with accelerated motion such as that at Mach 1 for uniform motion ([9]). 



3.2 Perturbers with positive acceleration 



To an accelerating perturber of trajectory, ^(t) = ^ot + /2 with ^o > 0, corresponds the phase A/r — ^ — r/2 where 
^ = t + ^0 is the Mach number. It could therefore seem sufflcient to take the force expressions of perturbers initially at 
rest and substitute ^ for t. There is however a key difference between the two cases: the condition from the retarded sonic 
wavefront, r ^ t, refers to time regardless of the type of motion. It should therefore be replaced with r ^ ^ — ^o that 
we refer to as condition (iv). This difference makes the force dependent on the initial velocity in contrast to the expressions 
derived assuming a constant velocity ([9]). The physical reason for such a dependence was discussed in the Introduction: as the 
velocity varies, the far part of the wake that was launched in the medium corresponds to an earlier velocity. Combining the 
force from the far part of the wake to that near the perturber results in a force that has a memory of the earlier velocity. 

Inspection of the conditions (i,ii) of the previous section along with (iv) shows that the two cases ^o > 1 and ^o < 1 
may be considered separately. In the following we do not derive the force expressions for the counterparts of Cases I and IV 
as their space extensions vanish when ro tends to zero. 

For an accelerated perturber with an initial supersonic speed {^o > 1), the boundary conditions (i, ii, iv) away from the 
transient phase (Case I) are identical to those of Case V (|27|) except that in and r2, t is replaced with and the outer 
integration domain is r2 ^ r ^ — Provided that ./^ > 1 + [(1 — ^o]^ + 2ro]^^^ that we call Case VI, the friction force 
is given as: 



ro.^o(4 + ri - r|') + \J(o{rl + r{) + ro(ri - r{ - 4)].^ + \rl + r| + J(o{r{ 



r2 



' + irl 



+ 



2^ + ro 
2^3 



log 



2J( + ro 



r\J( ■ 



ro 



■ ro 



2^3 



log 



2^ — ro 


2 , 


./# — ^0 




r-ro 


J( + J^O 



(30) 



New terms depending the initial Mach number appear in the expression of the force. As ./#o is set to zero the force reduces to 
H28|) with t — This shows that including acceleration in the perturber's motion makes the medium store the information 
about the perturber's initial state which is in turn reflected in the reaction force felt by the perturber. Taking the limit ro <C 1, 
the force becomes: 



Dynamical friction for accelerated motion 7 



Three extra terms appear with respect to the uniform velocity expression (|9]). As ^ ^q, the force tends to that of 
Case V (|29|) of perturbers initially at rest as the density perturbation of a larger Mach number dominates the initial density 
perturbation. 

For an accelerated perturber with an initial subsonic speed, the boundary conditions (i,ii,iv) show that the force expression 
is given as; 

1- [(l-.^o)'-2ro]"'' l + (2r-o)'/' : I ^ hi, (32) 

l + (2ro)'/'s;.^^l+[(l-.^o)'-2ro]''": / = 7iii, (33) 

1+ [(l-^o)' + 2r-o]'''' / = /vi, (34) 

where in the expressions of lu and Jm, t is replaced with We recover the general statement we made about the force on 
subsonic motion in Case II. In particular, the friction force on an initially subsonic perturber depends on the initial velocity 
as soon as the velocity is supersonic and becomes larger than 2c — Vq. In Fig. (4) we show the friction force with ro = 10"* for 
an initial subsonic velocity with = 0.6 using equations (|32|) . (|33|l and (|34p . and an initial supersonic velocity with = 2 
using (l30l) . 

We remark that the dependence of dynamical friction on the initial velocity of the perturber is not an artifact of choosing 
to assemble the perturber at t = 0. This starting condition can be replaced with t — —To and To in turn can be made to tend 
to infinity thus approximating a steady state. Examination of the general force expression shows that To ^ oo is equivalent to 
the condition t = where the limits ^ 1 and ^ <C 1 are taken for accelerating and decelerating perturbers respectively. 
For positively accelerating perturbers (Case IV), the force is independent of the initial Mach number as 3> 1 In the 

next section, we show that for decelerating perturbers the dependence on remains as the velocity decays to zero. 



3.3 Decelerating perturbers 



For decelerating perturbers < 0), we seek the dynamical friction force from t — until the velocity vanishes at t — \^o\. 
The choice of a negative Mach number instead of a negative acceleration makes the derivations easier but results in a positive 
friction force (opposed to the direction of motion). The phase A/r = ^ — r/2 is now negative and the boundary conditions 
that determine the integration domains of ((Sjl are (i), (vi) and (v) r ^ /(r, ro). The latter condition shows that Tmax = 1. 
Three cases need to be considered: —1 ^o, —1 — (1 — 2ro)^''^ ^ ^ —1 and ^ —1 — (1 — 2ro)^'''^. For subsonic motion 
( — 1 ^ ^o) and away from the transient regime (i.e. —1 + [(1 + ./#o)^ + 2ro]^/^ ^ the integration domains are given as: 



^ r ^ r: 



with 
with 



= -1. 



(35) 
(36) 



As for Case II (equations 1131 and 1 14p . the latter domain (|36p does not contribute to the force. The corresponding integrand is 
proportional to 1 — sign(l + — r/2) and the sign argument is positive for r < 2(^ + 1). The contribution of domain (|35p is 
identical to that of Case II and / — lu (|15p with t replaced by Dynamical friction on a decelerating subsonic perturber is 
therefore independent of its initial velocity. This confirms the general statement we made about subsonic motion in an infinite 
medium in Case II. 



For supersonic perturbers with an initial Mach number —1 — (1 — 2ro 



.1/2 



-1, we need to consider three cases 



depending on the value of the Mach number 
^ ^ —1 + (2ro)^^'^ whose boundaries are: 



The first which we call Case VII corresponds to — 1 — [(1 + ./£o)^ — 2ro]^''^ ^ 



< r < r. 



s; r s; 



«o, 



with Tniin = /(r, ^,ro), where r{ = . 

with Tmin = —1. 



■+l + [(.^ + l)' + 2r( 



11/2 



(37) 
(38) 



Domain (I37p contributes a term similar to (llip with Vi and r( substituted for t and r( respectively. Domain (|38p contributes 
a term similar to (|20p with r( and ^ — substituted for and r2 respectively. The force is then given as: 



/vii 



2^ + ro , 



2^ + ro 



r(.y£ + ro 



, 2.# - ro , 



r^ + rf)].^ 


' + 








+ (rf - rj)^^ 


o^^{^ + 


Jit, 


) 








2Jt — ro 


+ 


2 , 


J( — ^0 






r\Ji + n, 









The second case which we call Case VIII concerns — 1 + (2ro)^^^ ^ 
integration domains are: 



ri ^ r ^5 r, 
J 



2, with Tmin = /(r, .-#,ro). 



with 



-1, 



(39) 

-1 + [(1 + - 2ro]^''^ The corresponding 

(40) 
(41) 
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with 
with 



= -1. 



where 



' + 1 + [(.^ + 1)' - 2ro] 



1/2 



(42) 
(43) 



The integration radii satisfy < rj ^ rg ^ 2(./# + 1) ^ r{ . In contrast to Case III, this ordering shows that the second 
domain (|4H) does not contribute a finite term to the force. Domain (|40p contributes a term equal to In, and domain (|42|l a 
term similar to (|20p where rg and r{ are substituted for f and r| respectively. The last domain (143 [I contributes a term similar 
to that of psp in Case VII. Summing up all terms, we find: 



^o''o(r( -r{ +r{ - rl- 4) + [.Mo(rl + r{ + r| ~ r^) + ro(4 - rj; + - + r| 



[r\ + r-f + - + ./^o(r-f + r{ 



■ n 



^ + {r( + r{- ri - rj) 



+ 



2^ + ro , 



2./^ + ro 



+ ro 



+ 



2^ 



■ ro 



2^3 



log 



2^ 



- ro 



rg^ 



ro 



+ 



■log 



2^ 



+ 



■ ro 



2^3 



log 



+ ro 



■ ro 



(44) 



The third case concerns —1 + [(1 + + 2ro]^''^ ^ ^ whose boundaries turn out to be those of subsonic motion 

^0 > —1. The corresponding force is therefore given by (|15p where t is replaced Fig. (5, left panel) shows the force on 
a perturber with an initial velocity = —1.7 using equations (|39p . (144 |l and (I15[l alongside the uniform motion friction 
force (|9}. In the supersonic regime the uniform motion force overestimates the drag by about 10% at its maximum. In the 
subsonic regime the force on accelerated motion departs significantly from the uniform velocity force before joining it at 
~ 2 + ^0 = 0.3 to assume the expression (llSp . In particular, when the velocity vanishes so does the force. 
Lastly, for a supersonic perturber with an initial velocity ^ —1 — (1 — 2ro)^''^, we find that: 



.1/2 



1 - [(1 + .^o) - 2ro]'/^ «; < -1 + (2ro) 

-1 + (2ro)^''^ s; < 



I = Ivii, 
I = Iviii- 



(45) 
(46) 



Fig. (5, right panel) shows the force on a perturber with an initial velocity ^/£o — —4 using equations (|45|l and (|46p alongside 
the uniform velocity friction force ((9)l . In the supersonic regime the uniform motion force overestimates the drag by about 20% 
at its maximum. In the subsonic regime the force on accelerated motion departs significantly from that on uniform motion 
and is finite at ^ = 0. Taking the limit where ro is small, we find /vii(ro <S 1) = — /vi(ro ^ 1) and: 



i(ro < 1) 



4(1 + ./^o + ■ 



+ 



log 



+ .4 



2(1+. 



(47) 



At ^ = 0, /vni(ro -C 1) = 1 — 4/^o in contrast to the uniform motion force which vanishes with the perturber's velocity. 
The knowledge of the history of the perturber's motion is therefore required to evaluate dynamical friction. Simply applying 
uniform motion dynamical friction force ([9]) may yield incorrect results regarding the perturber's evolution in the ambient 
medium as evidenced in Fig. (6) where we plot the relative error between the accelerated motion and uniform motion forces 
for decelerating supersonic perturbers. A similar conclusion was previously reached by Sanchez-Salcedo and Brandenburg 
(1999) who numerically tested the accuracy of the uniform motion friction formula (|9} using the full nonlinear fluid equations. 
Following the evolution of a decelerating initial perturber with = 1.5 up to ^ = 0.6, they found that equation ([9| 
underestimated the friction force as ^ became smaller than unity (see their Fig. 5) which corresponds to our Cases VII and 
VIII. The authors attributed this discrepancy to the 'memory of the wake'. 



4 CONCLUDING REMARKS 

In this paper, we studied the influence of acceleration on dynamical friction in a gaseous medium. Our main conclusions 
are stated in the Abstract. This work also constitutes an illustration of the use of the general expression of dynamical 
friction derived in Appendix A. In particular, applying that formulation to a perturber with accelerated motion allowed us to 
illustrate how dynamical friction is not a local force and may depend the perturber's initial state. The meaning of the initial 
state depends on the problem's conflguration. For instance, it may correspond to the time the perturber enters the gaseous 
medium or to the initial time of a numerical simulation of dynamical friction. We note that non-locality is not peculiar to 
rectilinear motion and should be true of the dynamical friction associated with rotational (circular or eccentric) motion. 

The main technical result is the set of expressions h to /viii that describe dynamical friction as a function of the initial 
and current velocities as the perturber is accelerated (positively or negatively) in the ambient medium. Our derivations have 
assumed a constant acceleration because the calculations are simpler to carry out in this case. Dynamical friction for constantly 
accelerating perturbers applies readily to the motion of young stars with asymmetric jet systems. We believe that the general 
features we found will hold for a general type of acceleration as long as the instant velocity change remains small. Including 
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the deceleration generated by dynamical friction instead of a constant negative acceleration is likely to further reduce the 
friction force because the acceleration's amplitude increases with time. 

The dependence of dynamical friction on the porturbcr's size, ro, is most important for supersonic motion. In the force 
expressions, the dimensionless radius, ro, stands in for the ratio roA/c^ of dimensional parameters. The reaction force therefore 
depends explicitly on acceleration even in the limit ro ^ /A as in Case VI. Dynamical friction on subsonic motion is 
independent of m when ro -C /A. Using the force expression requires the knowledge of the porturbcr's acceleration and the 
medium's sound speed. In Fig. (3) to (5), the value ro = 10"* was used so as to make accelerated motion friction closer to 
uniform motion friction (see Fig. 1). 

The assumption of constant velocity in the derivation of dynamical friction is also made for coUisionless star systems. 
Although they have different magnitudes, the uniform motion friction forces in a gaseous medium and a coUisionless star 
system behave similarly (as \F\ oc for supersonic motion and \F\ oc V for subsonic motion; see Ostriker 1999, Sanchez- 
Salcedo and Brandenburg 2001). This similarity suggests that the conclusions we found are likely to hold for coUisionless 
systems. 
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APPENDIX A: GENERAL EXPRESSION OF DYNAMICAL FRICTION 

In this appendix, we derive a general expression of dynamical friction on a perturber of gravitational potential 4>p with a 
general trajectory ^(t) that travels inside an inviscid homogeneous gaseous medium. The perturber is set in motion at time 
t = 0. The medium's response to the perturber's presence is described by the standard fluid equations: 

dtp + V ■{pv) = 0, dtv + {vV)v = --Vp-je{t)V(l>p, (Al) 

P 

where v{x,t) and p{x,t) are the perturbed velocity and density fields, p is the pressure, and Jff{t) is the Heaviside function. 

Prior to the porturbcr's motion, the gaseous medium is assumed to be at rest [i.e.v(x,t) — 0], homogeneous of density po, 
pressure po and isothermal with a sound speed c = (po/po)^^^. The linearization of medium's equation in the vicinity of the 
rest state yields: 

dtg + cV- 1^ = 0, ^dtv + VQ=-'^^V4>v, (A2) 

where we used p = po{^ + o) and v = cu. The combined equations in turn yield the forced sound wave equation satisfied by 
the density perturbation g as: 

V'Q-^dig=-^V''cp,. (A3) 

The solution of this equation is known in electromagnetism theory as the Lienard-Wiechert potential (Jackson 1998, Chapter 
14) and represents the electric potential of a moving electric charge distribution. It is written in integral form using the 
retarded Green function as: 



In reaction to the perturber's excitation of the wake, the gaseous medium exerts a force on the perturber given as: 



1 f 

— OO 

reaction to the perturber's 

F = [ pogJ^m4>.d'x = I A'x<^'x'du ^'^-V^^' ' f (")] ^^^^t" ' -^(") (A5) 

■/av(.) 47rc2 Jg^^^j |x-x'| 

The symbol dV{x) denotes the boundaries of the spacial integration domain defined by the medium's physical boundary 
£is well as the lower distance cutoff in the perturber's vicinity such as its size or its corresponding accretion radius. The 
boundaries dV{x) are used to truncate the perturber's potential at small and large distances to avoid divergence. Applying 
the variable changes, X = x — £,{t) and 1^ = x — x' and integrating the force over u yields the general expression of dynamical 
friction as: 

^ ^ ^^mpoj d^X^^ P,(X-Y + A)VMX), (A6) 

where pp is the perturber's density distribution, pp = V'^(f>p/4-!TG, and A = ^{i) — ${t — \Y\/c). For a given perturbing 
potential, dynamical friction therefore depends only on the phase A as well as the boundary conditions specified by dV 
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and ^ ct. The dependence of A on the whole trajectory ^{t) shows that dynamical friction is manifestly non-local. 
Deriving dynamical friction through equation (|A6|I is in contrast with the previous two-step methods of estimating dynamical 
friction for uniform motion by first calculating explicitly HA4|) then applying an analytical or numerical calculation to the 
first integral in equation (IA5|I to derive the force (Ostriker 1999, Kim and Kim 2007, Kim et al 2008). As a result such 
methods obscure the true dependence of the force on the perturber's motion and the boundary conditions. For a point-like 
perturber, 4>p{x) — ~GAd/\x\ and Pp{x) = M5{x), where M is the perturber's mass. The forced sound wave equation, density 
enhancement and friction force are reduced to: 

1 1 It GM , 

V'e ~ ^ dh = -^V^0 = -^^^ - ^Wl-^W' (A7) 
,ix,t) = GM r--5[u-t+\x-iiu)\/c]J^iu) 



■je{t){GMfpo 



rsmedrded<p -2 — — — , (A9) 

dVly+i{t-r/c)],r!ict 'V ~ ' 



where y = x — ^{t — r / c) . The condition on the time of perturbation (t ^ 0) represented by the Heaviside function in equation 
(|A3|) and the expression (|A5|) now implies that the spacial domain that enters the expression of the force is equal to the size of 
the wavefront, ct. However, the spacial domain refers to the relative position with respect to the perturber not at time t but 
at the retarded time t — r/c. This explains why the force (|A9P describes subsonic as well as supersonic motion that perturbs 
the medium outside the sonic wavefront of radius ct. In the force expression, the boundaries dV[y + ^{t — r/c)] are treated as 
truncation boundaries. In the next Appendix, the force expression is applied to a point-like perturber in uniform rectilinear 
motion to recover the corresponding known results. 



APPENDIX B: APPLICATION TO UNIFORM RECTILINEAR MOTION 



We validate the expressions (|A6|I and (|A9|) by deriving dynamical friction for a perturber set in uniform rectilinear motion in 
an infinite homogeneous medium. The perturber's trajectory is written as ^(t) = Vte^ where is the unit vector along z, 
the phase A = ./Mre.^^ where .M is the Mach number. The friction force is equally along and its expression is simplified as 
follows: 



F = 



F 



.^{t){GMfpQ 



r sin 6Ar A9 dip 



r cos 6^ — A 



dV\y+(,(t-r / c)],r^ct 



2tv Jf{t){GMfpo 

V2 



dr 



8V[y+i(t-r/c)].r!ict 



|j.2 + A2 - 2rAcos6i|3/2 ' 
\-J(t 



(1 + _ 2.^T-)l/2 



(Bl) 
(B2) 



where in the last expression we performed the integrals over and r = cos Q. The extremal values Tmin and Tmax are determined 
from the boundary conditions. Calling ro the perturber's size, the force is truncated at small distances by |x — > ro. 
Expressing these conditions in terms of the variables r and r that appear in equation (IB2|) shows that we need to find the 
domains that satisfy: (i) —1 ^ r ^ 1, (ii) ro 5S r(l -f — 2^tY^'^ , and (iii) r ^ ct. The first two conditions yield the 
following integration domains: 



\l + yf 
ro 



<r < 



< r 



ro 



with Tmin = — 1, 
with Tmin = — 1, 



1 + .^2 _ (^^^^^2 

2^ 



(B3) 
(B4) 



The inner domain (|B3P corresponds to the vicinity of the perturber's retarded position. The radius ro/\l + ^\ defines a region 
around the perturber that does not contribute to dynamical friction as the first two conditions are not satisfied. The third 
condition depends on the distance ct travelled by the retarded wavefront. There are two force regimes: the local perturbation 
regime {ct ^ ro/\l — where the perturbation remains in the perturber's vicinity and the distant perturbation regime 
(ro/jl — ^\ ^ ct ) where the perturbation has expanded beyond the perturber's vicinity. Calling / the integral that appears 
in (IB2I) . the integration for the distant perturbation regime extends over the two domains HB3[) and (|B4p with r ^ ct and / 
is given as: 



I = 



/ = 



1 
2^ 



' log 



dr — [1 — sign(l 
] [1 + sign(l - J( 



dr 
r 



[1 — sign(l — 



log 



(1 - .£)ct 



ro 



(B5) 
(B6) 
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This expression agrees with those derived by Ostriker (1999) for subsonic and supersonic perturbers moving in an infinite 
gaseous medium using explicit expressions of the density perturbation. In particular, the force, F oc 7/./#^, is time-independent 
and behaves as for small Mach numbers. We note the interesting result that only the vicinity of a subsonic perturber's 
retarded position (IB3|l contributes to dynamical friction. We show in section 3 that this property implies that the force on 
subsonic perturbers with a general type of motion (not necessarily rectilinear or uniform) is similar to (IB6|l provided that the 
perturber's motion history does not involve supersonic episodes and that the perturber's acceleration \^{t)\ is much smaller 
than (? Itq. 

For the local perturbation regime (ct ^5 ro/ll — >^|), the force is zero before t — ro/c\l + ^\ and only the domain (|B3|I 
is relevant to the integral / which is written as: 

/ --0 2ro I ^ + ^2 ) ^ 

In Fig. (7) , we show the friction force as a function of t for subsonic and supersonic perturbers as t increases smoothly from 
to a value in the distant perturbation regime. The time dependence of the force agrees with that obtained from the numerical 
simulation of the nonlinear fluid equations (Sanchez-Salcedo and Brandenburg 1999, Figure 2). However, the attempted fit of 
the force in the local perturbation regime from numerical simulations (Sanchez-Salcedo and Brandenburg 2001, Equation 22) 
does not reproduce the correct dependence of (|B7p on ro and time. 



dr 



log 



ct(l -h. 



(1--^- 
2ct 2ro 



(B7) 
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Figure Bl. Evolution of the density enhancement q from a perturber initially at rest. The level curves are determined for values of 
Q = log(l + c'^q/GMA). Top left panel: subsonic motion, (1 ^ ^ ^ 3, = 0.1). The density increases nearer the perturber. Note how 
the level curves do not possess front- back symmetry with respect to the perturber. Top right panel: motion at Mach 1. Middle left panel: 
supersonic motion inside the sonic wavefront. The solid level curves correspond to (0.5 £) ^ 2, A§ = 0.1) and (2 < ^ sC 3, = 0.2) 
near the perturber's Mach cone. The dashed level curves correspond to a shallow depression around § = 0.527. Middle right panel; the 
perturber reaches the sonic Shockwave. The solid level curves correspond to (0.5 ^ § ^2, = 0.1). The dashed level curves correspond 
to shallow depressions around g = 0.37 inside the sonic wavefront and g = 2.3 inside the Mach cone. Bottom left panel: supersonic motion 
at = 4. The solid level curves correspond to (0.6 ^ g ^ 1, A^ = 0.1). The dashed level curves correspond to a shallow depression 
around g = 0.22. Bottom right panel: supersonic motion at = 8. The solid level curves correspond to (0.2 ^ ^ ^ 0.6, A^ = 0.1). The 
dashed level curves correspond to a shallow depression around g = 0.08. For supersonic motion, larger values of g are confined to the 
sides of the Mach cone where the discontinuity with the unperturbed medium is largest. 
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Figure B2. Evolution of the density enhancement g from decelerating perturbers. The level curves are shown for a subsonic perturber 
with ^0 = —0.5 (left column) and a supersonic perturber with = —2 (right column). Note how front-back symmetry is lost 
everywhere for the subsonic perturber. This however does not significantly modify dynamical friction from its value obtained using a 
constant velocity. 
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Figure B3. Friction force on an accelerating perturber initially at rest for three values of ro = 10^^, 10~*, and 10~^. The solid curves 
correspond to the friction force that takes acceleration into account given by equations II12II . I|15| l. I|21|l . J25} and (|28j, and the dashed 
curves to the uniform motion friction force expression | [9l . 




Figure B4. Friction force on accelerating perturbers with an initial subsonic velocity (left panel, ^ = 0.6) and an initial supersonic 
velocity (right panel, .,#0 = 2). The solid curves correspond to the friction force that takes acceleration into account given by equations 
II3UI I. 1 1328 . fi'il and 1134 I I. and the dashed curves to the uniform motion friction force JSJ. The small cutoff radius is tq = 10~^. The dotted 
curve in the left panel is the friction force on a perturber initially at rest. 
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Figure B5. Friction force on decelerating perturbers with initial supersonic velocities |./#o| = 1.7 (left panel), |./#ol = 4 (right panel). 
The solid curves correspond to the friction force that takes acceleration into account given by equations l|15^ . I|39^ and ||44| | for |^o| = 1-7 
and by equations 1 145 I I and I I46II for |^ol = 4. The dashed curves correspond to the uniform motion friction force l[9}- The small distance 
cutoff is ro = 10~^. 
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Figure B6. Relative error, 100 X |1 — |^o/^vn|| from equations ^ and l|39| l. on the friction force when using the uniform motion 
expression instead of the accelerated motion expressions for decelerating perturbers with initial supersonic speeds. The small distance 
cutoff is ro = 10~^. The curves correspond to different initial velocities. 
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Figure B7. Friction force for uniform motion as a function of time. The force, F, is given by Imt . llB6t . and (|B7ll and scaled to 
1-k(GM')^ I (?\ time is scaled to ro/c. The perturber's motion is initiated at t = 0. The Mach number is indicated to the left of each 
curve. The time dependence is similar to that obtained from the simulation of the non-linear fluid equations in Figure 2 of Sanchez-Salcedo 
and Brandenburg (1999). 



